Superconducting films with antidot arrays - novel behavior of the critical current 
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Novel behavior of the critical current density jc of a regularly perforated superconducting film is 
found, as a function of applied magnetic field H . Previously pronounced peaks of jc at matching 
fields were always found to decrease with increasing H . Here we found a reversal of this behavior for 
particular geometrical parameters of the antidot lattice and/or temperature. This new phenomenon 
is due to a strong "caging" of interstitial vortices between the pinned ones. We show that this 
vortex-vortex interaction can be further tailored by an appropriate choice of the superconducting 
material, described by the Ginzburg-Landau parameter k. In effective type-I samples we predict 
that the peaks in jc{H) at the matching fields are transformed into a step-Uke behavior. 



I. INTRODUCTION 

For practical applications of superconducting (SC) ma- 
terials, the increase and, more generally, control of the 
critical current in SC samples are of great importance. 
In recent years much attention was given to the investi- 
gation of superconducting films patterned with a regular 
array of microholes (antidots), which have a profound in- 
fluence on both the critical current and the critical mag- 
netic fieldiiSiSi^. Due to the collective pinning to the 
regular antidot array, vortices are forced to form rigid 
lattices when their number "matches" integer and frac- 
tional multiples of the number of pinning sites at fields 
Hn = n^o/S, and Hp/g = ^^q/S (where n,p,q are in- 
tegers) respectively, where $o is the flux quantum, and 
S is the area of the unit cell of the antidot lattice. This 
locking between the pinning array and the vortex lat- 
tice is responsible for the reduced mobility of the vor- 
tices in applied drive and consequently the increased crit- 
ical current at integer and fractional matching fields, was 
confirmed both by experiments (imaging^, magnetization 
and transport measurements-'^'^) and molecular dynamics 
simulationai. 

However, regardless on the imposed pinning profile, the 
vortices at interstitial sites always have high mobilitji^, 
show different dynamic regimes from the pinned ones, 
and their appearance is followed by a sharp drop in the 
critical current^. In this respect, the maximal occupation 
number of the antidots (saturation number ris) becomes 
very important for any study of the critical current. In 
an early theoretical work, Mkrtchyan and Schmidt^ have 
shown that rig depends only on the size of the holes. How- 
ever, in the case of periodic pinning arrays this number 
depends also on the proximity of the holes, on tempera- 
ture and on the applied fieldii^iiSiii. 

Besides the pinning strength of the artificial lattice, 
vortex-vortex interactions are crucial for vortex dynam- 
ics. Most of the experiments on perforated superconduct- 
ing films are carried out in the effective type-II regime 
(k, = 2K'^/d > I/V2, with K being the Ginzburg- 
Landau (GL) parameter, and d being the thickness of 
the SC film scaled to the coherence length ^), where vor- 
tices act like charged point-particles, and their interac- 



tion with the periodic pinning potential can be described 
using molecular dynamics simulations (MDS)^. However, 
the overlap of vortex cores (with sizes ~ ^) , and the exact 
shape of the inter-vortex interaction (depending on the 
material properties refiected through k), which arc ne- 
glected in MDS, may significantly modify the equilibrium 
vortex structures and consequently the critical current. 

In the present letter we investigate the critical current 
jc of superconducting films with regular arrays of square 
antidots, taking into account all parameters relevant to 
the SC state, within the non-linear Ginzburg-Landau 
theory. This formalism allows us to analyze the jc de- 
pendence on the geometrical parameters of the sample, 
material (even type-I) and temperature, and compare our 
results with existing experiments. 



II. THEORETICAL APPROACH 

We consider a thin superconducting film (of thickness 
d) with a regular lattice of square holes (side a, period W) 
immersed in an insulating media with a perpendicular 
uniform applied field H (see the inset of Fig. To 
describe this system, we solved the nonlinear Ginzburg- 
Landau (GL) equations for the order parameter ip and 
the vector potential A (in dimensionless units, and with 
temperature T taken explicitly in units of the zero field 
critical temperature Tco, see Refji^ for more details): 
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The magnitude of the applied field H is determined by 
the number of flux quanta piercing through the sim- 
ulation region Ws x Wg- On the hole-edges we used 
boundary conditions corresponding to zero normal com- 
ponent of the superconducting current. Periodic bound- 
ary conditions''^ are used around the square simulation 
region: A{r + li) = A{r) + VCi(r) and ■0(r -I- U) = 
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FIG. 1: (Color online) The critical current density jc (a) (in units of jo = cHc2£,/4'jtX^) as a function of the applied magnetic 
field H (in units of the first matching field Hi) for different antidot-sizes and fixed period W = 8^, and the contour plot of the 
Cooper-pair density at H — H2 (b) and a.t H = Hz (c) for R = 2.6^. The insets show the schematic view of the sample and 
the antidot occupation number Uo as a function of the antidot-size a at different matching fields (i.e. for 2 — 5 vortices per unit 
ceU). 



-0 exp(27r«Ci(r^/(f>o)7 where li {i — x,y) are lattice vec- 
tors and (i is the gauge potential. We use the Landau 
gauge ^0 = HxCy for the external vector potential and 
Cx = HWs-y + Cx, Cy = Cy, with Cx, Cy, being con- 
stants. Generally speaking, depending on the geometry 
of the antidot lattice, one must minimize the energy with 
respect to the latter coefficients. The size of the supercell 
in our calculation is typically 4x4 unit cells (containing 
16 holes, i.e. Ws = AW). We solved the system of Eqs. 
(P)-© self-consistently using the numerical technique of 
Ref:^^. In order to calculate the critical current, first 
we determine the ground vortex-state for a given applied 
magnetic field after multiple starts from a randomly gen- 
erated initial Cooper-pair distribution. Then the applied 
current in the x-direction is simulated by adding a con- 
stant Acx to the vector potential of the applied external 
field^i. Note that the current jx in the sample result- 
ing from the applied Acx is obtained after integration 
of the x-component of the induced supercurrents (calcu- 
lated from Eq. over the y = const, cross-section 
of the sample. With increasing A^x^ the critical current 
is reached when a stationary solution for Eqs. (Q-© 
ceases to exist (i.e. vortices are driven in motion by the 
Lorentz force). 



III. INFLUENCE OF GEOMETRICAL 
PARAMETERS OF THE SAMPLE 

The enforced stability of the vortex lattices in period- 
ically perforated superconducting sample a^i^^ at integer 
and fractional matching fields leads to pronounced peaks 
in the critical curren t^i^i^ . However, the exact shape 
of these lattices, and consequently their stability when 



locked to the pinning arrays, strongly depend on the pa- 
rameters of the sample. For example, while small anti- 
dots can pin only one vortex, in larger holes multi-quanta 
vortices may become energetically preferable This re- 
duces the number of interstitial vortices, whose higher 
mobility affects strongly the critical current of the sam- 
ple. Therefore, in this section we investigate the critical 
current of our sample for different sizes of the antidots a 
and antidot lattice period W . 

Fig. n] shows the critical current density jc of the sam- 
ple as a function of the applied magnetic field H for three 
different values of the antidot size: a = 1.5^, a = 2.6^ and 
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FIG. 2: (Color online) The critical current density jc of the 
superconducting film at matching fields for different periods 
of the antidot lattice. The inset shows the antidot occupation 
number Uo as a function of W at different matching fields. 
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a = 4^. The lattice period is W" = 8^, the film thickness 
d = 0.1^ and GL parameter equals k = 1 (corresponding 
roughly to Pb, Nb, or Al films). Note that, following 
the suggestion of Wahl^^, we took into account the in- 
fiuence of the perforation on the effective GL parameter 
K* = 2K^/d{l - 2a^/W'^). As shown in Fig. □ in the ab- 
sence of applied field the samples with smaller antidots 
always have larger jc simply due to more superconduct- 
ing material, i.e. larger screening. However, for H ^ 0, 
the critical current depends on the vortex structure in 
the sample. For small antidot-size a = 1.5^, where only 
one vortex can be captured by each hole (see the inset of 
Fig. the jc{H) curve shows the expected maxima at 
integer matching fields Hi, H2, H3 and H4 and at some 
of the fractional matching fields. As observed before^, 
the peaks of the critical current at matching fields de- 
crease with increasing magnetic field. However, a novel 
phenomenon is found for 2.5 < a/^ < 2.8: the critical 
current at the 3'"'' matching field is larger than the one at 
the 2"'' matching field. Moreover, the critical currents for 
H2 < H < H3 are higher than those obtained for fields 
Hi<H < H2. 

The explanation for this counterintuitive feature of jc 
lies in the hole occupation number rio, i.e. the num- 
ber of vortices inside the hole, and, consequently, in the 
saturation number of the holes (n^ — Uq for larger 
fields). For holes with 2.5 < a < 2.8, Ug is equal to 1 
(rig « a/4^^). Therefore, a.t H — Hi, all vortices are 
captured by the antidots. Analogously, at _ff = H2, be- 
sides the pinned vortices, one vortex occupies each in- 
terstitial site. If the same analogy is followed further, 
one expects two vortices at each pinning site at the 3'"'* 
matching field. However, this depends also on the dis- 
tance between neighboring holes. Namely, the period 
W determines the distance between the interstitial vor- 
tices. In other words, for small W, this distance is small, 
vortex- vortex repulsion is large, and it may become ener- 
getically more favorable for the system to force one more 
fiux-line into the hole rather than have two vortices at 
interstitial sides (see Figs.nib,c)). Therefore, our results 
show that the occupation number Hq of the holes in a lat- 
tice depends not only on the size of the hole a, but also 
on the period W and the number of vortices per unit- 
cell of the hole-lattice. For the parameters given above, 
Uo = 2 at H = i/a, i.e. only one vortex is located at 
each interstitial site, as for the case of = H2. Note 
that this interstitial vortex interacts repulsively with the 
pinned ones. This interaction is roughly twice as large at 
H = H3 than at H = H2, and the interstitial vortex be- 
comes effectively "caged" . This "caging" effect has been 
found for 7^ <W < 8.3^ for the radius R = 2.6^. For 
W > 8.3^, Uo and Ug become independent of W and 
for W < 7£_, the larger suppression of superconductivity 
around the holes at H — H3 becomes more significant 
than "caging" , and jc reverses again in favor oi H = H2. 

With further increase of hole-size a, more vortices are 
captured by the holes. Consequently, in order to observe 
the "caging" effect, one needs to consider higher magnetic 



fields, i.e. for Uq = 4, at least H > H^ is needed. In any 
case, this effect can always be realized at a given magnetic 
field by an adequate choice of a and W. Fig. |21 shows 
the critical current density of our sample for different 
periods of the antidot lattice W at a = 2.6^. For clarity, 
we plotted only the critical current at the matching fields. 
At H = and H = Hi (when there are no interstitial 
vortices) the critical current is an increasing function of 
W, because of the stronger screening by larger quantities 
of the superconducting material. For higher magnetic 
fields interstitial vortices nucleate more easily in a sample 
with larger W, which leads to a reduction of jc- The 
dependence of Uq (and thus also saturation number n^) 
on W is shown in the inset of Fig. El Uo decreases from 
Uo = 5 to Uo — 1 at H = H^, when we increase W from 
45 to 10^. 



IV. TEMPERATURE DEPENDENCE OF THE 
CRITICAL CURRENT 

So far, we presented results of our calculations at fixed 
temperature, where all units were temperature depen- 
dent (e.g. distances were expressed in C(T')). In what 
follows, we consider a superconducting film with thick- 
ness d — 20nm, interhole distance W = 1/im, and antidot 
size a = 0.33/L(m. We choose the coherence length ^(0) ~ 
40nm and the penetration depth A(0) = 42nm, which 
are typical values for Pb films. We study the infiuence 
of temperature with the help of the right-side term of 
Eq. (P) , which actually describes the temperature depen- 
dence of the coherence length ^(T) = £,{0)/ ~ T/Tcq\ 
and penetration depth A(r) = A(0)/^|1 - r/Tco|. 

Fig. O shows the calculated critical current density of 
the sample as a function of the applied field normalized 
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FIG. 3: (Color online) The temperature dependence oi jc{H) 
for a superconducting film containing an antidot array, at 
temperatures T/T^o = 0.86 — 0.98. The period of the antidot- 
lattice is W — l/im, the antidot size a = 0.33/im, film thick- 
ness d = 20nm, and k equals 1.05. 
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to the first matching field at temperatures T/Tco =0.86, 
0.88, 0.9, 0.92, 0.94, 0.96 and 0.98. As discussed in 
the previous section, the jc{H) curve shows pronounced 
maxima at matching fields with a substantial drop af- 
ter the number of vortices per unit-cell exceeds the hole- 
saturation number (in our case, for H > Hi). 

As expected, decreasing temperature results in an in- 
crease of the critical current for given magnetic field. 
However, the qualitative behavior of the jc{H) char- 
acteristics changes. Although the W/a ratio remains 
the same, the size of vortices and the occupation num- 
ber of the holes may change with temperature (as ^(T) 
changes). Actually, we found the same vortex structure 
at all temperatures T > 0.8Tco, but jdHz) is found to 
be larger than jc{H2) only in the temperatures range 
0.8Tco < T < 0.93Tco. At higher temperatures vor- 
tices become larger and suppression of superconductivity 
around the holes "masks" the critical current enhance- 
ment by the "caging" effect. This decrease of Cooper-pair 
density around the holes is also responsible for the de- 
creased jc{Hi)/ jc{0) ratio with increasing temperature. 
At temperatures lower than O.STco, the peaks of jc again 
decrease with magnetic field, since W/^ significantly in- 
creases and its influence on Uq diminishes (see previous 
section). 

Interestingly enough, this jc{Hn+i) vs. jc{Hn) revers- 
ing behavior was recently found experimentally, but not 
noticed. In Fig. 6 of RefA a clear enhancement of the 
critical current at H/Hi — 3 was found which is larger 
than the one at H/H2 as is similar to the behavior found 
in Fig. 13 for T/Tco < 0.94. A quantitative compari- 
son between theory and experiment is difficult because 
of the different determination of jc- In our calculations 
we assume normal state as soon as vortices are set in mo- 
tion, whereas in transport measurements a certain value 
of threshold voltage determines the critical current. 



squares) at T = 0.97Tco. For d = lOnm, the sample is 
still in the type-II regime (k, = 2.887), and the criti- 
cal current shows a peak-like behavior at the matching 
fields. The drop in the critical current for H > Hi is 
caused by the interstitial vortices which is larger with in- 
creasing K*. Note that the "caging" effect is also present 
for these values of the parameters. When the film thick- 
ness is increased {d = 50nm, k» — 0.577) the critical 
current density is higher for H < H2. For H < Hi, this 
increase is achieved due to a stronger Meissner effect. At 
fields Hi < H < H2 the increase of jc is more appar- 
ent, as all vortices are captured by the holes (except for 
d = lOnm, where Uq — 1). As soon as interstitial vortices 
appear in the sample {H > H2), jc becomes even smaller 
than the one for d = lOnm. This inversion of jc clearly 
demonstrates the higher mobility of interstitial vortices 
in type-I superconductors. 

Looking at the jc{H) curve as a whole for d = 50nm, 
we observed a pronounced step-like behavior. Note that 
in type-I samples the matching between the number of 
flux-lines and the number of antidots does not lead to 
a peak- like increase of the critical current. Namely, re- 
gardless of their number, additional flux lines are doubly 
pinned by the attractive hole-potential and the attractive 
interaction with previously pinned vortices. The "step" 
in jc{H) occurs only when the number of vortices per 
hole Uo changes (more vortices and consequently a larger 
suppression of if) around the hole), or when interstitial 
vortices appear. The effect of an increase of Uo diminishes 
as d increases (i.e. decreases and screening increases), 
or when temperature is lowered (as discussed before). 
This tendency is illustrated by solid and open dots in 
Fig. 21 and ultimately leads to a two-step jc{H) curve, 
with larger critical current for H < H^, and smaller jc 
for H > Hn, where n — Hq. 



V. EFFECTIVE TYPE-I VS. TYPE-II 
BEHAVIOR OF THE CRITICAL CURRENT 

In the previous sections we have shown that higher mo- 
bility of interstitial vortices leads to a dramatic decrease 
of the critical current. Therefore, keeping the interstitial 
sites of the superconductor vortex-free is essential for an 
improvement of jc. In this respect, we consider effec- 
tively type-I superconductors, where (i) the screening of 
the magnetic field is enhanced {i.e. vortices will be com- 
pressed in the holes), and (ii) the interaction of vortices 
becomes attractive (depends on lnK»). As an example 
we considered a sample, with antidot size a = O.bfim, 
lattice period W = 1.5/im, ^0 = 40nm and Ao — lOnm. 
We fine-tuned the effective vortex-vortex interaction by 
varying the thickness of the sample, i.e. changing from 
type-II to type-I behavior with increasing d. Fig. 2]shows 
the critical current of the sample for four values of the 
film thickness: d = lOnm (solid dots), d = 50nm (open 
dots), d = lOOnm (solid squares) and d — 150nm (open 
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FIG. 4: (Color online) The critical current density jc as a 
function of the applied magnetic field for different thickness 
of the sample. Remaining parameters are listed in the figure. 
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VI. CONCLUSIONS 

We studied the critical current of a superconducting 
film containing an array of antidots in a uniform mag- 
netic field, as a function of all relevant parameters. We 
found that the well-known jc enhancement by artificial 
vortex pinning strongly depends on the antidot occupa- 
tion number rio- The latter is determined not only by 
the size of the antidots as commonly believed, but also 
by their spacing and the applied field. As a consequence, 
when the parameter conditions are met, the critical cur- 
rent becomes larger at higher matching fields, contrary 
to the conventional behavior. Such a feature is a result 
of the "caging" of interstitial vortices between the larger 
number of pinned ones. This effect is strongly influenced 
by temperature, which agrees with a recent experiment.* 



Additionally, the interactions in this system can be tai- 
lored by the effective Ginzburg-Landau parameter . In 
effectively type-I samples, vortices attract, jc increases 
and peaks at matching fields diminish due to the enforced 
pinning at all fields. As a result, a novel step-like behav- 
ior of the critical current is found, with a sharp drop at 
higher fields, when interstitial vortices appear. 
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